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TECHNICAL FIELD 

The present invention relates to a method for measuring 
5 the Polarisation Mode Dispersion or PMD of an optical 
fibre, to the system able to allow its application and to 
the device able to perform the PMD measurements according 
to said method. 
BACKGROUND ART 

10 As is well known, polarisation mode dispersion is a 

phenomenon of distortion of the optical signal that 
propagates in single mode optical fibres, i.e. optical 
fibres that carry signals in a fundamental propagation mode 
constituted by a pair of degenerate modes with orthogonal 

15 polarisation. 

This phenomenon is, as is well known, linked to the dual 
degeneration of the fundamental mode of the optical fibres. 
Construction imperfections in real fibres (core 
ellipticity, external stresses, etc.) entail different 

2 0 group velocities for the two degenerate modes of the 
optical signal transmitted at one end of the fibre and 
hence distortions of the signal received by a receiving 
device (receiver) located at a second end of the fibre. 
Therefore, it is important to develop accurate techniques 

2 5 for measuring PMD to characterise, for example, 
telecommunication systems that use optical fibres to carry 
signals, in particular in the case of high speed or high 
capacity transmission systems, for instance at bit rates up 
to 10 Gbit/sec or higher. 

30 The PMD of the fibres is measured, for instance, by 

connecting a segment of optical fibre (fibre) 11 to be 
characterised to appropriate measuring instruments or . 
devices 10 (Fig. 1) . 



CONFIRMATION COPY 



wo 2004/005877 



PCT/EP2003/007048 



2 

The paper "Polarization Mode Dispersion of Short and Long 
Single-Mode Fibers" by N. Gisin et al, published on the 
"Journal of Lightwave Tech." Vol.9 on 7 July 1991 
discloses, for instance, an instrument 10 for measuring the 
5 PMD, which comprises an optical source 21 (Fig. 1, Fig. 2), 
able to generate optical signals to be applied to the fibre 
111 a Michelson interferometer (interferometer) 25 
connected to the second end of the fibre 11, able to apply 
Differential Group Delays or DGD to the received optical 
10 signal and, at the output of the interferometer 25, a 
measuring device 27 able to convert the optical signals 
into electrical signals and to process, on the basis of 
program modules stored therein, the electrical signals 
received, 

15 According to the prior art, the electrical signals 

obtained downstream of the interferometer 25 are 
graphically represented by the measuring device 27 by means 
of an interf erogram having in its x-axis time offset values 
introduced by the interferometer and in its y-axis values 

20 of intensity of the photo-current generated. 

In particular, according to the mentioned paper, said 
interf erogram is to be interpreted as a Gaussian 
distribution of PMD of the fibre, so that the PMD itself 
can be determined according to the standard deviation of 

25 the Gaussian curve that best interpolates the 
interf erogram. 

This type of interpretation of the electrical signals 
derives from the fact that to the fibre 11 is applied a 
propagation model called "waveplate model" , which consists 
30 of considering the real optical fibre as a great number of 
birefringent plates positioned in cascade with respect to 
each other and having their main axes oriented randomly. 
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This model is used, as is well known, to assess signal 
time-of -flight , in the hypothesis that at each interface 
between a waveplate and the next one the optical signal is 
coupled randomly on the main axis of the downstream 
5 waveplate. 

Based on the known model and in the realistic hypothesis 
that the length of the fibre is far greater than the 
coupling length between the polarisations, the prior art 
concludes that the measured graphic representation 

10 corresponds to a Gaussian curve with a standard deviation 
equal to twice the standard deviation of the optical signal 
due to PMD and hence indicative of the PMD of the fibre. 

However, it is believed that the measurements obtained 
taking the known model as a reference do not precisely 

15 correspond to real PMD values. 

In essence, it is believed that the conclusions of the 
prior art, which are based on the identification of the 
distribution of the times-of -flight of the optical signal 
with the distribution of the electrical field are wrong. 

20 It is also believed that the PMD measurements obtained 
according to the known model are consequently only 
approximately correct. 
DISCLOSXIRE OF THE INVENTION 

The technical problem the present invention aims to 

25 tackle is to determine a new propagation model of the fibre 
that will lead, in the first place, to a new method for 
measuring the fibre PMD that is more realistic under all 
computational conditions . 

The present invention therefore aims to solve the above 

30 technical problem by means of a new model for interpreting 
the interferogram and consequently a new method, system and 
device for determining the PMD of the fibre, as claimed. 
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In particular, the new model which is based on the 
theory of propagation of the optical signal as can be 
determined according to the so-called Principal States of 
Polarisation (PSPs) of the fibre, leads to a new 
5 interpretation of the interf erogram, characterised by the 
identification in said interf erogram of a central peak and 
two side lobes. 

Thanks to this new model, it is possible to determine 
the PMD of the fibre on the basis of the information 
10 content associated with at least one of the side lobes 
identified in the interf erogram , 
BRIEF DESCRIPTION OF THE DRAWINGS 

This and other characteristics of the present invention 
shall become readily apparent from the following 
15 description of a preferred embodiment, provided purely by 
way of non limiting example with the aid of the 
accompanying drawings, in which: 

Fig. 1 shows a general measurement layout that uses the 
interf erometric method (IM) ; 
20 Fig. 2 shows a block diagram of the system of Fig. 1; 

Fig. 3 shows a general diagram of the model of the optical 
source constituting the basis of the method according to 
the invention; and 

Fig. 4 shows an example of the results obtained with the 
25 method according to the invention. 

DESCRIPTION OF A PREFERRED EMBODIMENT 

With reference to Fig. 2, the instrument or system for 

measuring the PMD able to use the interf erometric method 

(measuring system or apparatus IM) 10 comprises, in 
30 particular, an interferometer 25 and a measuring device 27 

connected at the output of the interferometer 25 itself and 

able to control its operation, in known fashion, by means 

of a control connection 23. 
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The measuring device 27 comprises, for instance, a 
control unit (CPU) 30, a read-only memory unit (ROM) 31 and 
a random access memory unit (RAM) 33, all known and 
connected to each other in known fashion by means of a 
5 data, control and command channel (BUS) . 

The measuring device 27 further comprises a display device 
(display) 35, known in itself, and an opto-electronic 
module 37 associated to an amplifier 38, respectively able 
to convert the optical signals coming from the 
10 interferometer 25 into electrical signals and to amplify 
them; the display 35 and the amplifier 38 are connected, in 
known fashion, to the BUS. 

The ROM 31, according to the present invention, comprises 

program modules developed when designing the measuring 
15 system 10 and able to implement the method according to the 

invention, as will be described in detail farther on. 

Said program module, as will be readily apparent to a 

person versed in the art, can also be stored in the RAM 33 

by means of reading devices, not shown in the figure, able 
2 0 both to read from appropriate computer products the program 

modules according to the invention and to store said 

program modules in the RAM 33 itself. 

With reference to the measuring system 10, the method 

according to the present invention is based on the 
25 hypothesis of broad band radiation with spectral components 

that propagate according to the theory of the Principal 

States of Polarisation (PSPs) . 

In particular, with reference to the measuring system 

10, the method according to the invention is based on the 
30 following interpretative model and the corresponding 

logical -mathematical passages: 

(a) the spectrum of the (incoherent) optical source is 
approximated in segments by a family of rectangular 
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spectral intervals with extension Av (Fig. 3), each 
described by its amplitude, by its random 
polarisation state and by a phase factor, which is 
also stochastic; 
5 (b) the spectral width of the interval is specified 
relative to the band of the PSPs and to the 
bandwidth of the amplifier 38, in order to obtain a 
reasonably accurate approximation of the spectrum of 
the optical source 21, 
10 (c) a statistical description of optical propagation 
inside the fibre being measured is used, in terms of 
ensemble of spectral intervals defined above, 
ensemble that excites a statistical set of PSP, 

(d) the output spectrum of the signal is thereby 
15 obtained, which comes to depend on a set of 

statistical parameters: amplitudes of projection on 
the PSP and on the axes of the polarisation 
analyser, values of DGD and phase factors; 

(e) the interf erometric photo-current is expressed 
20 according to the power spectrum of the optical 

source 21 and of the density of probability of the 
DGD determined by the PMD of the fibre 11. 

The result obtained at passage (e) represents one of the 
characteristic elements of the present invention and is 
2 5 applied to implement the method according to the present 
invention . 

In particular, based on the proposed model, illustrated in 
Fig. 3, the spectrum S,-„(v) of a randomly polarised source 
is subdivided M-N rectangular spectral intervals S^{y) 
30 having width and amplitude U„ according to the 

following mathematical expression: 
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(1) 



Each replication translated in frequency of the fundamental 
interval has a polarisation state and a phase factor 
represented by statistical vectors p„ and by independent 
5 stochastic variables ip„ . 

The choice of spectral width 6^v is constrained by two 
extremes of a physical nature. The upper limit for Ai/ is 
the minimum band of the PSPs in the spectral domain of the 
optical source 21 (Fig. 2, Fig. 3), since the assumption is 

10 that every spectral interval propagates according to the 
model of the PSPs at the first order. This upper limit is, 
as is well known, in the order of the inverse of the 
average DGD for the fibre being measured and it falls 
within the GHz range even for values of DGD of tens of 

15 picoseconds. On the other hand, the lower limit for Au is 
imposed by the opto- electronic module 37 and by the 
amplifier 38 which, as is well known, has a limited band, 
as set out in detail hereafter and it can fall in the range 
of hundreds of kHz . 

20 A value of Av selected within the frequency domain 
identified by said upper and lower limit is to be 
considered suitable for an approximation in rectangular 
segments and sufficiently accurate for all LED optical 
sources commonly used in IM apparatuses (10) ; moreover, 

25 this conclusion confirms that the approximation proposed in 
the mathematical expression (1) has correct physical 
foundations . 

The set of spectral intervals { UnSe(v- nAv) } , after 
being launched in the fibre 11, excites a set of PSP pairs, 
30 each characterised by the value zlr„ of DGD. At the input 
end of the fibre, the elementary spectral intervals are 
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projected on the slow and fast PSP, respectively with 
amplitudes ^ and ^^-a„ ; in which a„ represents the power 
coupling coefficient between the n^^ interval and the axis 
associated with the corresponding fast PSP. 

In similar fashion, at the output of the fibre 11, the 
coefficients of projection of the output PSPs on the axis 
of polarisation of the interferometer are indicated as ^ 



and yj^-Pn . The spectrum of the signal after the latter 
selection in polarisation is then expressed in the form: 

in which, for the sake of convenience, the two new 
statistical variables a„=^a„fi„ and b„ = ^J^^^^ajly^^ have been 
introduced. 

The average value of the DGD of the PSP was neglected, 
15 referring the quantities Ar„ to the origin of the delays 

set at the interferometer. 

Since in the time domain, as is well known, the 

interf erometric photo-current I is proportional to the 

auto-correlation function of the output signal, in the 
20 frequency domain the quantity to be computed is the power 

spectrum : 

= I Z ""n^nS, - nA v)e-'''^'^''e'^^ + f^b,U,S, {v - nA vy'^'-^'e'^' ] • (c.c.) 

(3) 

relationship in which 'c.C indicates the complex 
2 5 conjugated expression of the term in round parentheses. 

According to the present invention, it is hereby 
demonstrated that the interference between different 
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spectral intervals Sq(v -nAv)e^^" does not give rise to any 
significant contribution to the photo-current and therefore 

can be neglected. 

For this purposes, the expression of the photo-current 
5 generated by the opto -electronic module (detector) 37 and 
by the amplifier 38 shall now be discussed in detail, with 
particular reference to the band limitation of the detector 
37 and of the amplifier 38. 

The module of the field incident on the opto-electronic 
10 module 37 has, as is well known, the following expression: 

in which 7' = 2-(A, -L2) is the delay introduced by the offset 
c 

of the interferometer, Ep (t, T) is the electrical field 
incident on the opto-electronic module 37 and Ejn (t) is 
15 the electric field of the signal at the output of the 
fibre . 

It should be noted that all signals depend on two 
variables: the variable t is the time variable whereon the 
signals are measured, with a domain coinciding with the 

20 time interval of the measurement (typically a few seconds) . 
The variable T is, as stated, the delay introduced by the 
offsetting of the interferometer 25, with a domain that 
coincides with the DGD interval of interest, typically 
[0.1, 100) ps. The quadratic photo detection and electrical 

25 amplification process is characterised, neglecting the 
conversion constants of the physical magnitudes, by the 
following relationships between the signals: 

hi' J) = 4(t.T) ^^[£in(0 + £iA( + T)f = ^[EUO + EUt + r) + 2£,^(0£w(/ + T)] 
I(t, T) = (r. 7-) • h(t) = ^[EI (/) + £f^ (/ + r) + 2Ej^ (t)E,^ (/ + T)] ♦ h{i) 
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in which XR(t, T) is the current generated by the opto- 
electronic module 3 7 and sent to the amplifier 38, I(t, T) 
is the current at the output of the amplifier 38 and h(t) 
is the impulse response of the amplifier 38 with limited 
5 band, which naturally depends only on the variable t. 

Recalling that the signals Ejn(t) and Kji/^t+Tj are at 
optical frequencies and that no electrical amplifier is 
capable of following their variations, one obtains: 

10 in which the symbol < >i/b indicates the time average 
computed over an interval in the order of the reciprocal of 
the band of the 'amplifier 38. 

Since the continuous component does not carry the 
information and is suppressed by the measuring instruments 
15 of the PMD, it is best to consider only the alternate 
component of the previous signal : 

^c(^n«£/^(0^/w(' + n♦A•(0 
in which the spectrum of h' (t) coincides with that of h(t), 

20 except in the origin where it is zeroed. 

In the frequency domain, computing the Fourier transforms 
with respect to the variable t, one obtains: 

fAc(AT) oc ^'(/) (4 ) 

It is readily apparent that the spectral components of the 

25 convolution in parentheses that exceed the band of the 
amplifier 38, are cut by the amplifier 38 itself. 
It is clear that, since the signal Ejulf) has very broad 
band, typically some THz, to proceed with the analysis it 
is necessary to subdivide the spectrum into sectors: 

30 f^S,(f-n^r)e''^' . 

Replacing this expression in (4) , one obtains: 
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According to the present invention, it is now possible 
to demonstrate that, under typical conditions of 
measurement and with commonly used amplifiers, thanks to 
the ratio between the band B of the amplifier and the 
extension Af of the spectral sectors, the expression (5) 
can be reduced to the diagonal form. 

Considering, then, the convolution of two adjacent 
rectangular spectral sectors ri(f) and rzlf) of amplitude A 
and band Af\ 

Hereinafter the convolution of an elementary spectral 
sector with itself will also be considered: 

To evaluate the influence of the convolution of adjacent 
sectors on the photo -current, one can compute the ratio 
between the energies of the two preceding convolutions in 
the band of the amplifier 38: 

B a B 

SNR^^-i = ^-5 B ' 

Jlr,(/)*r2(/)]2rf/ \{A^ffdf \[A'f^W 

0 0 0 



SNR = - 
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3 


B 
0 



A'B^ 



Given that k^Af/B, one easily computes: 

SNR^Zk^ -3it + 1 . 

This value is further divided by two, taking into account 
that each sector has two adjacent sectors. The following 
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Table 1 shows the value of SNR/2 thus computed as a 
function of the Af/B ratio. 



Af/B 


SNR/2 (lin.) 


Sm/2 (dB) 


2 


3.5 


5.45 


3 


9.5 


9.79 


4 


18.5 


12 ,68 


5 


30.5 


14 . 78 


6 


45.5 


16.59 


10 


135 . 5 


21 .32 



Tab. 1 



It is readily apparent that for k=3 the value of SNR/2 
5 already approaches 10 dB and exceeds 22 dB for k=10. 

One can therefore conclude that, having appropriately 
chosen the ratio k^Af/B, the prevalent contribution to the 
convolution in (2) is given by the sum of the convolutions 
of the homologous terms, and hence the current photo- 
10 detected by the opto-electronic module 37 can be computed 
as : 

Using the result just obtained, only the diagonal terms of 
the equation (3) survive and the photo-current I{v) takes 
15 on the following form: 

(7) 

The equation (7) is, according to the present 
invention, the most general expression of the photo-current 
20 in the frequency domain that can be obtained without making 
hypotheses on the statistical nature of the PMD of the 
fibre being measured. Hereafter, the two cases of fibre 
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with deterministic PMD and with random PMD shall be 
examined separately. 

An optical fibre with a deterministic birefringence 
(for example, a polarisation-maintaining fibre) is 
5 characterised by constant DGD and PSP as a function of 
frequency. Consequently: Ar„ - At , a„ = a , b„=b are constant 
and the equation (7) is simplified in: 

i.e., 

10 /(v)cc(a^ +6^)|5Jv)|' ^ab[e'^'^' ^e^"^%„{vf . 

A Fourier ant i -transform directly yields the expression of 
the photo- current in the time domain: 

I{T) ai{a'^b')s, (T) + ab[s, (I - Ar) + 5, + Ar)] ( 8 ) 

in which Sc(D is the source auto-correlation function, and 

15 corresponds to the Fourier ant i- trans form of |Sf„(v)|^ . 

The photo-current has three components: the source auto- 
correlation function and its replications translated about 
the two symmetric points ±Ar . This well known result [1] is 
obtained with simple calculation within the proposed model. 

20 Therefore, thanks to said characteristics the side 

lobes of the interf erogram have the same shape as the 
central peak and are located in the two points ±Ar , 
symmetrical with respect to the central peak. 

The case of fibre with random birefringence is the more 

25 important one in practice and requires a more complicated 
mathematical treatment. The analysis is founded on the 
computation of the average value of the random process 
represented by the photo-current I(v) . 

This computation has dual importance. First of all, it 
30 leads to a simple and elegant formal result represented by 
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equation (9) . Secondly, from the point of view of 
application to PMD measurement, this mathematical operation 
corresponds to a computation of averages of different 
measured interf erograms, as will be shown hereafter, 
5 allowing an experimental verification of the proposed 
model. From equation (7) one computes the statistical 
average (.) of the photo - current : 

Taking into account the linearity of the statistical mean 
10 and the fact that the random variables Ar„ are independent 
from the random variables aj}„ and [al-^b]) (this is a direct 
consequence of the random nature of the birefringence of 
the fibre), one immediately obtains: 

15 Since the statistical average of the exponential of the 
random variables Ar„ is their characteristic function 
M ' previous equation can be rewritten in the form: 

20 Moreover, for a given fibre, the characteristic functions 

^^,^{y) and the statistical averages {ajb„) and {al^bl) do not 

depend on frequency; i.e. they do not depend on the sum 
index n. Therefore, the previous equation is simplified as: 

25 A Fourier ant i- transform provides the statistical average 
of photo -current in the time domain: 
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(9) 



in which Ar(^) is the probability density of Ar , and the 
asterisk indicates the convolution. 

In this case, too, the photo-current is given by three 
5 contributions. The first term of equation (9) is linked to 
the properties of coherence of the optical source, and it 
is represented by a peak located in the origin of the times 
axis. The second and third term correspond to two 
symmetrical lateral lobes and contain all the information 

10 required to measure the PMD. The expression of the average 
photo -current, in spite of the symmetry about a peak in 
7 = 0, is very different from a Gaussian, due to the 
presence of the side lobes constituted by the convolution 
of s^(T) with the probability density of the differential 

15 group delay /^rC^) • Since it is assumed that the coherence 
time of the source (i.e. the width of the central peak) is 
negligible relative to the average DGD of the fibre, the 
two side lobes represent approximate measurements of the 
probability density of the DGD of the fibre. The quality of 

2 0 this approximation depends, as stated, on the coherence 
time of the source used, and can be improved by means of an 
appropriate de- convolution of the int erf erome trie diagram 
measured with the auto-correlation function of the source 
itself 3c(T) (this function can be measured separately by 

25 directly connecting the optical source to the 
interferometer) . Lastly, it can be stated that equation (9) 
is a generalisation of equation (8) if one only considers 
that the probability density of the DGD of a fibre with 
systematic birefringence is a Dirac function 5, 

30 Therefore, thanks to this characteristic the 

measurement of PMD with interf erometric technique is 



wo 2004/005877 



PCT/EP2003/007048 



16 

conducted not with a Gaussian interpolation of the 
interf erogram, but rather by analysing at least one of the 
side lobes of the interf erogram, considered, according to 
the present invention, as an approximate measurement of the 
5 probability density of the DGD associated with the PMD of 
the optical fibre 11. 

The propagation model according to the present 
invention was tested by means of simultaneous PMD 
measurements using both the interf erometric method, 

10 described in the present embodiment, and the method, also 
known, of JME or Jones Matrix Eigenanalysis . 

The tested fibre (50 km of coil wrapped single mode 
conventional fibre) was placed in a thermal chamber at 
temperatures varying from about 2«'C to about 60^C at 5®C 

15 steps. After each variation, the fibre was allowed to 
stabilise before conducting the interf erometric and 
polarimetric measurements. Thereby, over several days, 15 
interf erograms and 15 DGD measurements with JME method 
(each constituted by 91 points) were acquired. The data 

20 were then processed with the following procedure: 

• the mean of the interf erograms was computed after they 
were normalised; 

• said mean was compared with two symmetrical replications 
of the DGD probability density computed from the JME 

25 data. 

Note that the first processing operation is suggested 
directly by the interpretation of the interf erograms as 
probability density and corresponds to the computation of 
the mean photo-current as shown in equation (9) . 
30 From Fig. 4, which shows the results of the measurements, 
the good agreement between the results is readily apparent. 

In particular, from the results of the measurements, it 
is readily apparent thait the central peak of the 
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interferogram must be neglected because it depends only on 
the characteristics of coherence of the source and not on 
the PMD of the fibre. 

The result obtained, according to the present 
5 invention, attributes a new physical meaning to the 
interferometric curves, particularly useful for measurement 
activities . 

Obvious modifications or variations are possible to the 
above description, in the dimensions, shapes, materials, 
10 components, circuit elements, connections and contacts, as 
in the details of the circuitry and of the illustrated 
construction and of the method to be implemented, without 
thereby departing from the spirit of the invention as 
specified by the claims that follow. 

15 



